Spin-triplet pairing instability of the spinon Fermi surface in a f/(l) spin liquid 
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two dimensional spin liquid state. This phase is described by a two-dimensional spinon Fermi sea coupled to a 
t/(l) gauge field. We study Kohn-Luttinger-like pairing instabilities of the spinon Fermi surface due to singular 
interaction processes with twice-the-Fermi-momentum transfer. We find that under certain circumstances the 
pairing instability occurs in odd-orbital-angular-momentum/spin-triplet channels. Implications to experiments 
are discussed. 
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Understanding possible phases of matter in strongly inter- 
acting electron systems is one of the central issues in con- 
densed matter physics. A prime example of such systems is 
the Mott insulator. Among a plethora of possible insulating 
phases, perhaps the most interesting states are spin liquids 
with no long range order. The possibility of such states of 
matter was first discussed by Pomeranchuk back in 1941 |1], 
who conjectured that in insulators, the elementary excitations 
may be charge-neutral fermions (spinons). More precise and 
modern form of such proposals had to wait until Anderson 
proposed the spin liquid state as the key paradigm in high-r^ 
cuprates physics f^. 

While the relevance of spin liquid phases to the high- 
Tc problem is still under debate, the search for spin liq- 
uids has continued in other classes of strongly interacting 
electron systems. One of the most notable examples is a 
series of recent experiments on the organic compound k- 
(ET)2Cu2(CN)3 il li]. Here the electrons are half-filled on 
the triangular lattice and the relative interaction strength is 
controlled by applying a hydrostatic pressure. The insulating 
phase of this system exhibits finite uniform susceptibility and 
finite specific heat coefficient at low temperatures, i.e. metal- 
like behaviors. It was proposed that the U(l) spin liquid phase 
with a spinon Fermi surface is a viable explanation of the 
experimental results Si. Variational calculations on the 
Heisenberg-ring-exchange model are consistent with this ex- 
pectation [5J. 

The spinon Fermi sea, however, is strongly coupled to a 
U{1) gauge field. On the lattice, the U(l) gauge field is com- 
pact and in principle one should worry about possible confine- 
ment effect in two dimensions due to monopole events |7]. 
It has been argued that the monopole events are suppressed 
due to the coupling to the gapless spinon degrees of freedom 
0, Nevertheless, the gauge field still gives rise to sin- 
gular renormalization of various physical quantities. For ex- 
ample, the sgecific heat coefficient should diverge as r '^-^ 
12 , Tsl [14I1 . This behavior has not been seen and in- 
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stead there exist abrupt changes in the susceptibility as well 
as in the specific heat around a few Kelvin [3, 4]. This may 
suggest a thermodynamic transition or possibly a crossover. It 
is therefore interesting to consider possible instabilities of the 
corresponding spinon Fermi sea state. 



The instability of the spinon Fermi surface may also be a 
very useful avenue to study the emergence of other possible 
phases. Notice that in Fermi liquid theory, all broken symme- 
try states of ordinary metals can be understood as an instabil- 
ity of the Fermi surface. In the same spirit, the spinon Fermi 
surface state may be regarded as a mother state of various pos- 
sible phases of Mott insulators. 

In this paper, we study the Kohn-Luttinger-like pairing in- 
stabilities 1 15] of the spinon Fermi surface. Due to the cou- 
pling to the gapless gauge field, the spinon interaction vertex 
for the momentum transfer of twice the Fermi momentum di- 
verges as a power law [JJJ . This singularity renormalizes the 
effective interaction in the Cooper channel, which becomes 
attractive for odd orbital angular momenta. It is found that 
under certain conditions the system is unstable to the spinon 
paired state in odd orbital angular momentum (spin-triplet) 
channels. When the pairing is present, the ground state would 
be a Z2 spin liquid. 

Some remarks on the relation to the experimental findings 
are in order, a) As mentioned earlier, the singular temperature 
dependence expected from the gauge interaction has not been 
seen at low temperatures. This is consistent with the fact that 
the spinon pairing gaps out the U(l) gauge field below a pair- 
ing temperature so that various singularities associated with 
the t/(l) gauge field fluctuations do not show up at low tem- 
peratures, b) The specific heat does not depend on the applied 
field up to 8T 1 16]. This is consistent with the spin-triplet 
pairing that is not affected by the Zeeman effect, c) When 
the system becomes superconducting at higher pressure, the 
Knight shift does not change across the superconducting tran- 
sition 117]. This can be explained if the resulting supercon- 
ducting state is related to the spinon pairing state and hence a 
spin-triplet superconductor. 

Previously, Lee et al. lIlSll proposed a different mechanism 
for spinon pairing. It is based on the "Amperian pairing" that 
arises due to the attractive current-current interaction between 
spinons moving in the same direction. Such an interaction 
leads to Cooper pairs with a finite center of mass momentum, 
resulting in a translational-symmetry broken state. The main 
difference between our proposal and that of Ref. [6] is that 
our spin-triplet state occurs in the normal BCS channel and 
corresponds to a uniform ground state. 
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We start with the following Hamiltonian, which describes 
a system of fermionic spinons minimally coupled to a U(l) 
gauge field S E [H [B [11 Hi : 



r/;'(r)e(-/V-a)/,(r) 



(1) 



where e(-/V - a) is obtained by replacing p by (-iV - a) 
in the spinon dispersion e(p). Here and / are the spinon 
creation and annihilation operators, s - 1,2, ...A^ is the 
"spin" index generalized to components, and a is the fluc- 
tuating gauge field. The bare gauge field Hamiltonian can 

be taken as Ti, = 5 j d^r l^gle^(r) + 4 [V x a(r)]^| , where 

a and e = -/^ are canonically conjugate quantum operators, 
and go is a bare gauge coupling constant. However, the ex- 
act form and parameters of the bare gauge field theory are 
unimportant, as the coupling to spinons generates more rele- 
vant terms. This is pictorially described in Fig. 1.1a, where 
the dynamics of the field is determined by the particle-hole 
excitations of spinons via the RPA renormalization. The re- 
sulting effective propagator of the transverse gauge field takes 
the following form 



-iyoj/q+Xq'^' 



(2) 



where P^'^pi^) - ^afi ~ laqplq^, and y and x are con- 
stants, which determine the Landau damping and the Lan- 
dau diamagnetic susceptibility of the fermions, respectively. 
In the large-A^ limit, one can develop a self-consistent 
Eliashberg-type ifioil theory by using the electronic self- 
energy (Fig. 1.1b), S(e) = is\(i)Ql e\^^^ (where cdq ~ Ef is 
a constant). This leads to a non-Fermi liquid behavior. The 
self-consistency of the theory implies that (i) The lines in the 
bosonic and fermionic self-energies in Figs. 1.1a and 1.1b can 
be either thin (bare propagator) or thick (renormalized prop- 
agator), leading to the same result. Further renormalization 
does not change the effective dynamics of the fermions and 
gauge bosons, (ii) Once the above renormalization is taken 
into account, the vertex corrections are small (see Fig. 1.2). 
The latter statement is indeed true in the large-A^ limit, but 
only if the momentum transfer is not close to 2pf. As pointed 
out by Altshuler et al. lil2tl . the vertex diverges logarithmically 
if the momentum transfer is exactly equal to twice the Fermi 
momentum. The leading logarithms can be summed within 
the standard parquet technique leading to a power-law diver- 
gence of the vertex 



[\aj/Ef\ + c(q/2pF-lf^Y' 



(3) 



where To is the value of the vertex far from the 2/:if -anomaly 
and c and cr are some A^-dependent constants. In Ref. I1I2I1 . 
it was shown that within the large-A^ treatment, cr oc 1/A^. 
The extrapolation of the results in the large-A^ limit to N - 2 
leads to cr a; 0.36 for the circular Fermi surface. It can be 



shown that cr increases as the curvature of the Fermi surface 
becomes larger ^20^. On the other hand, the small -A^ approach 
in Ref. [12] gives cr ^ 0.52. Clearly, these results provide only 
an estimate of the exponent cr as there is no truly controlled 
method of treating strong gauge fluctuations. 
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FIG. 1: (la): The RPA approximation for the gauge-field propagator, 
(lb): The leading contribution to the fermionic self-energy. (2): The 
renormalization of the spinon vertex by gauge interactions. The cor- 
responding diagrams diverge logarithmically if the momentum trans- 
fer is equal to two Fermi momenta, q = 2pf. (3): The Bethe-Salpeter 
equation for the irreducible part of the scattering amplitude in the 
Cooper channel. The wavy line includes all interactions, including 
the magnetic interactions mediated by the gauge field. The shaded 
vertices represent processes described by Fig. 1.2. 

Here we assume that the structure of the theory found at 
large A^ is preserved in the physical limit of N - 2. This im- 
plies that the structures of all power laws in the bosonic and 
fermionic propagators are the same but the parameters of the 
effective theory can be some non-universal numbers. One of 
these parameters is the exponent cr in Eq. This param- 
eter also enters the temperature dependence of the suscepti- 
bility and therefore is a physically observable quantity. Thus 
instead of describing the effective theory in terms of the un- 
physical number of flavors A^, it is perhaps more reasonable to 
use the parameter cr as the variable which "controls" the ef- 
fective theory. This is the point of view we take in this paper 

We note that the power law divergence ([3]) can be viewed 
as a strong renormalization of the Kohn anomaly in the non- 
Fermi liquid phase. A natural question is whether there are 
instabilities related to this anomaly. Ref. j2l|l studied pos- 
sible instabilities in the particle-hole channel and concluded 
that a strong enough short-range interaction may lead to a 
density-wave transition in the system. Below we address the 
issue of possible instabilities in the Cooper channel. The the- 
oretical motivation comes from the following observation: In 
a usual Fermi liquid, the effective interaction with momen- 
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turn 2pp, is not divergent, but non-analytic. This weak non- 
analyticity is not benign and leads to a number of observable 
effects such as "non-Fermi-liquid" temperature corrections to 
thermodynamics and transport coefficients l:22i I23I |24 1 . This 
anomaly can also give rise to superconducting or pairing in- 
stabilities even in the presence of repulsive interactions. The 
latter phenomenon has been known as the Kohn-Luttinger ef- 
fect IT5, 25, 3. In a usual Fermi liquid, this effect arises 
because the effective interaction in real space acquires an os- 
cillatory part (Friedel oscillations), which is a reflection of the 
Kohn anomaly. Thus there appear attractive regions and elec- 
trons can pair up due to this attraction. They however, must be 
far apart from each other which implies a finite orbital angu- 
lar momentum of a pair In a gauge-fermion system, the 2pf 
anomaly is much stronger than in the ordinary Fermi liquid, 
thus it is natural to consider the possibiUty of spinon pairing 
of Kohn-Luttinger- type. 

Fig. 1.3 shows pictorially the Bethe-Salpeter equation for 
the effective interaction in the Cooper channel. Below we will 
concentrate on the case when the momentum transfer is of 
order 2pf. The wavy line in Fig. 1.3 is assumed to be a com- 
bination of the effective magnetic interaction and short-range 
repulsive interactions which may include the screened "elec- 
tric" forces mediated by the gauge field. From now on, we will 
denote the corresponding interaction as U{2pf). It is impor- 
tant to include the singular 2pf vertices in the Bethe-Salpeter 
equation. The resulting equation has the following form: 

Tcis,e';p,p') = U(s-e',p-p')rl (s-s',p-p') 



f/(e-77,p-k)r|,/e-77,p-k) 



X g,{m^,(-k)T, (77, s'; k, p') , (4) 

As usual in the Cooper problem, we assume that the integral 
on the right-hand side is determined by momenta at the Fermi 
surface and small frequencies ^ 0. Therefore, the mo- 
mentum dependence of all functions in Eq. (|4| reduces to the 
dependence on a single angular variable (i.e., the angle be- 
tween two of the following vectors p, p', and k). We there- 
fore can simplify the Bethe-Sapleter equation by performing a 
standard decomposition into harmonics corresponding to dif- 
ferent orbital angular momenta, i.e. fi - f {(j)) COS {l<p)d<p In, 
where /(0) is an arbitrary function of the angle. In these no- 
tations, we can write Eq. (|4|i as 

Tc{s,e'-l)^U{2pF)[Tl,^\{e-s\l) 

-U(2pF)T 2 [r^^J (s - ij, DC(T])r, (7?, , (5) 

n 

Here C(j]) is the Cooperon and we assumed that the interaction 
(without the vertex corrections) has no singularity atq - 2pf. 
Pairing instability arises in principle only if there is an effec- 
tive attraction in one of the channels labelled by the orbital 
angular momentum I. In the large-A^ limit, the parameter cr 
may be small and by performing the Fourier transform of the 



square of the 2pf -vertex, we obtain (cr < 1/6) 
K.] - -3a + 1,1 -3a -I) 



(-1) 



(6) 

where B{p, q) - Y{p + q)l \r{p)T{q)\ is the beta-function and 
the last estimate on the right-hand side corresponds to the limit 
Z » 1. Therefore, for small cr < 1/6, the effective interaction 
is attractive for large even orbital angular momenta. However, 
in the physical limit of N - 2, there is no reason for the pa- 
rameter cr to be small. If o" > 1/6, the corresponding Fourier 
transform of the static double-vertex diverges, but it is cut off 
by the frequency and we get (cr > 1/6) 



r , T 2rlB(l/3,2o-- 1/3) , 

3,,./3 (-1)' 
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In this regime, the effective interaction is attractive for all odd 

orbital angualr momenta, / - 1,3,5, 

We now address the possibility of a pairing instability due 
to the effective attraction. In the usual Fermi liquid an at- 
tractive interaction automatically implies a pairing instabil- 
ity, which is a consequence of the Cooper logarithmic di- 
vergence. It is not so in the fermion-gauge system because 
the Cooperon divergence is much weaker here than in the 
Fermi liquid [we recall that ^;,(k) - (/t/Iwq/'?!'^'' - ^p)^]: 

<^('?) = / (0^.(k)^-.(-k) = vn/(\?j\^/'ajl'\ where v is the 
density of states. 

To find an instability, we need to consider the following 
BCS self-consistency equation 

TU(2pF) 2 [r^^J (e - 77, C(77)<D,(77) = <D,(e), (8) 
1 

where I is the orbital angular momentum and O/ is related to 
the corresponding pairing amplitude. The existence of a non- 
trivial solution of the corresponding eigenvalue problem im- 
plies a divergence of the resolvent of the integral (if T = 0) or 
finite difference (if T + 0) equation (|5]i; this can be interpreted 
as a pairing instability (see also Ref. [26]). To unequivocally 
establish the existence of a pairing instability for given values 
of the bare interactions, we look for an infrared divergence in 
Eq. ([8]|. For this, we assume that the eigenvector is a weakly 
dependent function of the frequency. With this Ansatz, we get 



l+(- 



Ef 

-livKU{2pF)Ef-^l^ ^ ■ 



ds 



,|2o-+l/3 



O, = 0, 



where = (2/3)rf,B(l/3,2o--l/3)(£'/r/wo) ' is just a dimen- 
sionless constant, v - m/n, and we assumed that the Fermi 
energy serves as a high-energy cut-off. Clearly a non-trivial 
solution always exists if / is an odd number and cr > 1/3. 
The critical value cr = 1/3 gives a logarithmic divergence 
and the transition temperature T, ~ Ef exp {-l/[vKU{2pf)]}, 
while for larger values we get a power-law singularity and 

T, ~ £f [v7ff/(2pf)]l/(2tr-2/3)^ 
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We note that the transition temperature depends on the or- 
bital angular momentum very weakly (apart from the odd- 
even dependence); an explicit dependence on the value of I 
appears only in the subleading orders and enters via the over- 
all cut-ofF constant. Due to this weak dependence of the tran- 
sition temperature on I, the actual pairing symmetry is non- 
universal and depends on the structure of the effective inter- 
action far from 2pf and/or details of the fermion dispersion. 
What our argument above explicitly shows is that in the the- 
ory with cr > 1/3, the system is unstable to the pairing at some 
odd orbital angular momentum (or momenta). Thus the paired 
state is a Z2 spin liquid ll27ll . 

If the pairing occurs in a single odd angular momentum 
channel, then at some finite temperature of order T, we ex- 
pect a crossover from the spinon-Fermi surface phase to a 
pseudogap-like phase as the temperature decreases. On the 
other hand, the susceptibility and specific heat coefficients 
seem finite at low temperatures below the "kinks" at a few 
Kelvin. There may be two possible ways to explain this: 
a) The low temperature phase is not a homogeneous phase. 
Various extrinsic effects (such as tiny amount of disorder, 
structural inhomogeneity, and warping of the spinon Fermi 
surface) can lead to the suppression of the local transition tem- 
perature and the local pairing amplitude in certain regions of 
real space. The resulting state would consist of the gapless 
spinon puddles (where the pairing amplitude vanishes) in the 
sea of a Z2 spin liquid. This may be consistent with the ob- 
servation that two different NMR relaxation times have been 
observed L3j,|4J. The finite susceptibility and specific heat co- 
efficient at low temperatures (below the "kink" temperature) 
may be due to the gapless regions, b) Another possibility is 
that the pairing occurs in a large orbital angular momentum 
channel or the pairing amplitude is given by a superposition 
of contributions from multiple angular momentum channels. 
Then a large number of nodal excitations may exist and prac- 
tically it may look like there is an extended gapless region in 
momentum space. Which scenario is realized in real materi- 
als may depend on non-universal physics and it requires more 
detailed analysis of the band dispersion and sub-dominant in- 
teractions in the Cooper channel. This is beyond the scope of 
this paper, but it may be an excellent topic of future study. 

In conclusion, we demonstrated that there exist Kohn- 
Luttinger-like pairing instabilities in the t/(l) spin liquid with 
a spinon Fermi surface. This pairing occurs in odd-angular- 
momentum/spin-triplet channels. As explained in the intro- 
duction, this is consistent with several experimental findings. 
On the other hand, the previous proposal of the "Amperian 
pairing" |@] (while it explains various aspects of the same 
experiments) suggests a spin-singlet pairing and a spatially 
modulated pairing amplitude with a well defined wavevector 
This corresponds to an incommensurate version of the valence 
bond solid phase and would induce lattice distortion. Thus 
the difference between two proposals may be tested by X-ray 
scattering. 

Finally we emphasize that the microscopic scenario pre- 
sented here may be relevant not only to the gauge-fermion 



system of spinons but also to other fermionic systems with 
singular interactions, which include electrons near an Ising 
ferromagnetic instability or near a nematic ordering quantum 
critical point [12811 and the vortex metal phase proposed in 

Refs. ilHIiiF 
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